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We report a reversible and hysteretic change in the topography measured with a scanning tun¬ 
neling microscope near a single Fe dopant in a GaAs surface when a small positive bias voltage is 
applied. First-principles calculations of the formation energy of a single Fe atom embedded in GaAs 
as a function of displacement from the substitutional site support the interpretation of a reversible 
lattice displacement of the Fe dopant. Our calculations indicate a second stable configuration for 
the Fe dopant within the lattice, characterized by a displacement along the [111] axis, accompanied 
by a change in atomic configuration symmetry about the Fe from four-fold to six-fold symmetry. 

The resulting atomic configurations are then used within a tight-binding calculation to determine 
the effect of a Fe position shift on the topography. These results expand the range of demonstrated 
local configurational changes induced electronically for dopants, and thus may be of use for sensitive 
control of dopant properties and dopant-dopant interactions. 


I. INTRODUCTION 

The manipulation of single substitutional impurities in 
a host semiconductor has become increasingly important 
in the development of nanoscale semiconductor devices, 
permitting a wide range of tunable phenomena. Recent 
advances in nanofabrication techniques allow the precise 
placement of an individual dopant atom in a host lat¬ 
tice, as well as the construction of semiconductor device 
sizes such that a single atomic impurity can significantly 
affect the electronic, magnetic, and optical properties of 
the hostP, or form the active component of the devicd^. 
Scanning tunneling microscopes (STM) allow remarkable 
precision and detail that is ideal for investigating the 
electronic structures of these solotronic devices, and have 
been successfully used to map their spatial structures and 
wave functions with atomic resolution 4 . Recently, elec¬ 
tronic control of the lattice position of individual dopants 
using a STM has been demonstrated, from the reversible 
displacement of a single Si dopant in the surface layer of 
GaAs from a substitutional to an interstitial site^El, with 
a corresponding change of the dopant charge state. This 
transition is similar to the naturally occurring formation 7 
of a DX~ center in AbGai-^As alloys with x > 0.22. 

Magnetic, transition metal (TM) dopants in III-V 
semiconductors are of considerable interest due to the 
internal spin degrees of freedom associated with the core 
d-electron states. Control of the core occupation of indi¬ 
vidual TM atoms would make it possible to manipulate 
the core spin and sensitive magnetic interactions of the 
dopants. In recent worlPthe manipulation of the charge 
state of single Fe dopants at the [110] surface of GaAs 
using cross-sectional scanning tunneling microscopy (X- 
STM) was demonstrated. Direct electrical manipulation 
of the core occupation of a single Fe dopant from the 
isoelectronic 3d 5 (Fe 3+ )° configuration to the ionized ac¬ 
ceptor state 3d 6 (Fe 2+ ) _ was achieved using tip-induced 


band bending (TIBB) by the STM, changing the core 
spin state from S=5/2 to S=2. Under the conditions of 
Ref. 0 for a small positive bias voltage of the STM tip 
(where we normally observe a contrast related to the 3d 5 
(Fe 3+ )° configuration, see Fig. [lja)) we frequently ex¬ 
perienced a sudden switch of the contrast. The contrast 
after switching is shown in Fig. 0b). The rather peculiar 
electronic contrast of the topography in the proximity of 
the dopant consists of a bright strongly anisotropic fea¬ 
ture superimposed on a dark symmetric background. 

After a brief period, typically of ^ 10 seconds, the 
contrast switches back to the original contrast (which 
we refer to as type I) shown in Fig. [lja). We followed 
the switching behavior of single Fe atoms by repeatedly 
scanning over the same scan line in the [001] that is in¬ 
tersecting with the contrast of the Fe impurity. In Fig. [2] 
several time-traces are shown that consist of repeated 
scan lines taken in the [001] direction which are plotted 
next to each other which demonstrate the time depen¬ 
dence of the switching process. In all panels the switching 
appears to occur instantaneously, and the switch to the 
type II contrast shown in Fig.JTJb) lasts for about 10 sec¬ 
onds. The dark background observed for the type II state 
we interpret as due to negative charge that is localized 
at the scale of single atom, which affects the tunneling 
current!^. 

The shift seen in topography is not associated with a 
shift in the charge state alone. In Ref.[8]by modifying the 
occupation of the core d-level the charge state of the Fe 
atom was changed to 3d 6 (Fe 2+ ) _ . This negative charge 
state binds a hole in an effective-mass state with a bind¬ 
ing energy of 25 meV 10 . If the bright anisotropic contrast 
for Fe in the type II state were due to a hole bound by 
25 meV then a much larger scale and triangular shaped 
contrast should be visible, such as that associated with 
the 25 meV bound state of Zn 1 ^^, which is not the case. 
The bright anisotropic contrast suggests a much deeper 
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FIG. 1: STM topography images of the same Fe atom that 
is located about 3 ML below the (110) surface of GaAs be¬ 
fore and after switching. In (a) the expected (type I) con¬ 
trast for Fe in the isoelectronic 3 d 5 (Fe 3+ )° configuration is 
shown whereas in (b) type II contrast after the switching is 
shown. Both type I and type II topography were observed 
at the same small positive bias voltages. The type I topog¬ 
raphy is perfectly symmetric with respect to the [001] axis 
and highly symmetric with respect to the [110] axis. The 
type II topography consists of a bright anisotropic feature on 
a dark background. The dark background is attributed to the 
Coulomb field of negative charge located near the center of 
the Fe impurity in the type II state. Image size 7.5 nm by 8.7 
nm 


binding energy and therefore a relaxation of the lattice 
position of the Fe dopant from substitutional to intersti¬ 
tial lattice position may be responsible for this observa¬ 
tion. Such a switching scenario is similar to that which 
is seen in the case of Si in the surface layer of GaAs, 
in which the configuration of the Si atom changes from 
substitutional to interstitial with an associated change of 
charge stated. 

The existence of a second possible configuration would 
require two minima in the potential energy landscape 
in configuration space separated by a potential barrier. 
For the case of Si in GaAs, the configuration associated 



time 


FIG. 2: Time dependence of the switching behavior of indi¬ 
vidual Fe atoms. The time traces are obtained by repeatedly 
scanning over the same line intersecting with the topographic 
contrast of a Fe atom. In each panel the tip was repeatedly 
scanning along the same line in the [001] direction but the 
line is located at a different position in the [1-10] direction 
with respect to the position of the Fe atom. The length of 
the time traces is 120 seconds and the length of the scan line 
in the [001] direction is 9.5 nm. 


with the second minimum is described by a change of the 
charge state from the neutral sp 3 (Si 2+ )° state to the neg¬ 
atively charged sp 2 (Si + ) _ state. In addition, the charge 
density is strongly localized near the Si atom, causing a 
structural relaxation corresponding to an approximate 1 
A displacement along the [111] axis (along the axis of the 
Si-As bond) and away from the 4-fold symmetric substi¬ 
tutional position, resulting in a broken bond 7 . At the 
GaAs [110] surface, the formation energy of the intersti¬ 
tial configuration is negative and stable, while in bulk 
GaAs the formation energy is positive and metastable 6 . 
To evaluate the possibility of a reversible, electronically- 
induced shift in the lattice position of the Fe dopant we 
use first-principles calculations to investigate the elec¬ 
tronic structures for both cases of substitutional Si and 
Fe dopants in bulk GaAs and their associated formation 
energy profiles in configuration space. 


II. THEORETICAL & COMPUTATIONAL 
METHODS 


Following the previous work on the energetics of the 
DX~ center by Chadi and Chang^, the DX~ center is 
formed when an impurity in the 4-fold symmetric substi¬ 
tutional site, denoted by d, which gives up one electron 
to the conduction band, 
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FIG. 3: Energies of Si (a) and Fe (b) in configuration space embedded in bulk GaAs for the neutral and negative charge 
states. For convenience, the zero of energy is taken to be the total energy of the neutral charge state with the dopant atom in 
the 4-fold symmetric substitutional configuration, and the displacements along the [111] axis are relative to the substitutional 
configuration. 


(a) (b) 


FIG. 4: (a) 4-fold symmetric substitutional lattice position of 
Fe in GaAs corresponding with type I Fe. (b) Interstitial lat¬ 
tice position of negatively charged Fe in GaAs corresponding 
with type II Fe. 


d°^d + + e~, (1) 

where e is a free electron in the conduction band and the 
superscripts signify the charge state. The free electron 
then recombines with a different neutral impurity to form 
the DX~ center, 

<f + e~->DX~. (2) 

where DX is the broken-bond configuration. Thus, the 
complete reaction necessary for the formation of the 


DX center is given by the equation 

2 d°^d + + DX~. (3) 

The formation energy for the case of a negatively charged 
structure must then take into account the addition of the 
free electron from the conduction band. This is included 
by amending the calculated total energies of the structure 
by the difference of the energy of the conduction band 
minimum relative to the Fermi energy. 

In this work, we employ ground-state density func¬ 
tional theory (DFT) calculations using the linearized 
augmented plane wave (LAPW) method as implemented 
by the WIEN2k code 13 . The Perdew-Burke-Ernzerhof 
(PBE) generalized gradient approximation (GGA) ex¬ 
change correlation functional was used. For each self- 
consistent calculation a 7x7x7 k-point grid was used and 
the wave functions were expanded in a plane wave basis 
set with a cutoff energy of 10 Ry. The evaluated struc¬ 
tures consist of a 2x2x2 64-atom supercell of GaAs with 
a single dopant atom substituting in a Ga cation posi¬ 
tion. For the neutral and singly negatively charged cases, 
the dopant atom was displaced in increments of approx¬ 
imately 20 pm along the [111] axis away from the neigh¬ 
boring As atom to simulate the broken As bond, and the 
structure optimized at each step with the dopant position 
held fixed until the forces acting on each atom are less 
than 2 meV/A. For other charge states of Fe, the atom 
was displaced in increments of 40 pm along the [111] axis 
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with the structure optimized at each step until the forces 
acting on each atom are less than 5 meV/A. For calcula¬ 
tions performed with charged supercells, the final ener¬ 
gies were amended using the Makov-Payne correction to 
first ordei^ 4 . The Makov-Payne correction eliminates ar¬ 
tifacts which arise in ab initio calculations with charged 
supercells due to the periodic boundary conditions im¬ 
plemented in the calculation. This is performed by an 
addition to the final energies of the form 


where q is the additional charge, a is the lattice- 
dependent Madelung constant, and L is the linear size 
of the supercell. 

In order to calculate the effect of the shift from sub¬ 
stitutional to interstitial lattice position on the Fe wave- 
function, the ab initio results have been used as input 
for a tight-binding calculation. The tight-binding calcu¬ 
lation, which is an extension of the approach of Ref. [15] 
uses an sp 3 model for the GaAs host with the addition of 
d orbitals on the Fe dopant. In order to include the effect 
of the shifted position of the Fe, we use the generalization 
of Harrison’s d ~ 2 scaling law to calculate the new value 
of the hopping matrix element!^ for the shifted position. 
In the shifted case we amend the hopping between the 
Fe d orbitals and the three nearest As p-orbitals from 
their non-shifted values based on the ab initio calculated 
bond lengths. Hopping from the Fe to the original fourth 
nearest neighbor As atom is taken to be zero, simulating 
a broken bond. 


III. RESULTS AND DISCUSSION 

The energy profiles in configuration space for the neu¬ 
tral and negative charge states of both Si and Fe can be 
seen in Fig. [3j Initially the calculations described in the 
previous section were performed for the case of Si in bulk 
GaAs for both charge states. The location of the sec¬ 
ond minimum of the formation energy of the negatively 
charged Si in the DX~ center configuration was found 
to be at an approximate 1.37 A displacement from the 
4-fold symmetric substitutional site along the [111] axis, 
while the interatomic distance from the broken-bond As 
increases from 2.40 A to 3.47 A resulting in a net in¬ 
teratomic displacement of 1.07 A from the neighboring 
As atom, in good agreement with previous calculations^. 
The DX~ configuration of Si is one of 3-fold coordina¬ 
tion symmetry. The formation energy of the broken-bond 
configuration, Ef , is calculated in the same manner as 
found in the work of Yi et alP in which only energies be¬ 
tween same charge states are compared. The formation 
energy is calculated from the equation 

( 5 ) 


where d~ is the energy of the negatively charged dopant 
atom at the 4-fold symmetric Ga substitutional position 
and DX~ is the negatively charged dopant in the inter¬ 
stitial configuration of lowest total energy. Using this def¬ 
inition, the formation energy for Si in the DX~ configu¬ 
ration was found to be 353 meV, suggesting a metastable 
configuration as expected for the case of in bulk GaAs. 
Previous calculation^ have found that this formation en¬ 
ergy switches from positive to negative as the Si atom is 
brought to the surface layer of [110] GaAs. 

These calculations were repeated for the cases of neu¬ 
tral and negatively charged Fe in bulk GaAs. The result¬ 
ing structures corresponding with the energy minima for 
type I and type II are pictured in Fig. [4j For the case 
of negatively charged Fe, the location of the second min¬ 
imum of total energy in configuration space corresponds 
to an approximate displacement of the Fe atom by 1.76 A 
from the 4-fold symmetric substitutional position along 
the [111] axis, while the interatomic distance between the 
Fe and broken-bond As increases from 2.32 A to 3.87 A 
resulting in a net increase of 1.55 A. The associated for¬ 
mation energy of this configuration was found to be 1.20 
eV, implying that the interstitial configuration is only 
metastable in bulk GaAs, similar to the case of Si. The 
resemblance of the energy profiles for the Si and Fe cases 
in bulk GaAs is suggestive that the switching behavior 
of Si in bulk and near the surface may serve as an anal¬ 
ogous model for the behavior of Fe at the surface. This 
would imply that the formation energy of the negatively 
charged Fe in the broken bond configuration may become 
positive as approaches the surface layer, a behavior which 
would be supported by the experimental observation of 
a change of topography near the ionized Fe dopant near 
the surface layer of [110] GaAd^. 

One notable difference between the Fe and Si cases is 
the greater displacement of the Fe atom along the [111] 
axis by approximately 0.39 A, which places the atom into 
a configuration of 6-fold coordination symmetry as op¬ 
posed to the 3-fold coordination symmetry seen with the 
Si DX~ configuration. Another difference appears when 
comparing the neutral cases of Si and Fe. Whereas the 
neutral Si atom is completely unstable at all displace¬ 
ments away from the 4-fold symmetric substitutional con¬ 
figuration, the neutral Fe energy profile has a shallow sec¬ 
ond minimum, approximately 103 meV in depth, with a 
similar atomic displacement to the interstitial configura¬ 
tion for the negatively charged Fe. However, this second 
minimum for the neutral Fe has an associated formation 
energy of 1.79 eV, 590 meV greater than that of the neg¬ 
atively charged configuration. 

The doubly negative and singly positive charge states 
of Fe in bulk GaAs both exhibit a second minimum, how¬ 
ever they are found to be very shallow, similar to the case 
of neutral Fe, with depths of approximately 99 meV and 
59 meV for the singly positive and doubly negative charge 
states respectively. These energy minima correspond to 
displacements of the Fe along the [111] axis greater than 
2 A. The formation energies are found to be 1.82 eV and 


E f = E(DX~) - E(d~) 
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FIG. 5: Wavefunctions of the t 2 symmetry state of Fe three monolayers above the position of the Fe for (a) the substitutional 
position and (b)-(e) interstitial position with different symmetries relative to the (110) plane. The cartoons beneath the 
wavefunction plots depict the bond symmetry for each plot relative to the the (110) plane containing the Fe. 



FIG. 6: Topographic images of the same Fe atom in the type I 
and type II state. For both types the contrast is symmetric 
relative to the (110) plane containing the Fe. The strong con¬ 
trast features appearing along the [111] directions cross each 
other in the same point for both type I and type II. This cross¬ 
ing point coincides with the position of a Ga atom in the sur¬ 
face of GaAs. The dark background related to the Coulomb 
field of atomically localized negative charge is shifted by about 
a lattice constant along the [001] direction as indicated by the 
black circle and cross that are centered with respect to the 
dark background in the contrast. 


2.31 eV for the doubly negative and singly positive charge 
states respectively, with total energies greater than both 
the singly negative and neutral cases. These energy pro¬ 
files suggest that there are no viable configurations of 
the Fe dopant along the [111] axis which exhibit an as¬ 
sociated doubly negative or singly positive charge state, 
leaving the (Fe 2+ ) - acceptor state as the only metastable 
charge state for the interstitial Fe configuration. 


The tight-binding wavefunction calculations further 
support the origin of type II Fe as due to a shifted lattice 
position of the Fe dopant. In Fig[5^a) the bow tie shape of 
the type I (substitutional) Fe t 2 symmetry state is clear. 
In Fig [5jb)-(c) the bond between the Fe and a nearest 
neighbor As below/above the (110) plane containing the 
Fe has been broken and the Fe shifted towards its nearest 
neighbors in the [001] direction. In Fig|5jd)-(e) the bond 
between the Fe and one of its nearest neighbor As in the 
(110) plane has been broken and the Fe shifted closer to¬ 
wards its other nearest neighbors. In the tight-binding 
calculation the largest change to the bonding configura¬ 
tion between type I and type II Fe is the broken bond 
between the Fe and one of its nearest neighbor As, which 
does not destroy the overall bow tie shape of the Fe wave- 
function. 

Close inspection of the topographic contrast for the 
same Fe atom in either the type I or type II state shows 
that the bright and dark contrast features are symmetric 
relative to the (110) plane containing the Fe. The to¬ 
pographic images furthermore show that the strong con¬ 
trast appearing along the [110] directions, especially on 
the lefthand side of the images, cross in a position that 
coincides with a Ga atom at the surface. This can be 
understood if the Fe atom is located in an odd layer be¬ 
low the semiconductor/vacuum interface. We estimate 
that the Fe atom is in the 3rd layer counting from the 
semiconductor/vacuum interface. The symmetry center 
of the dark contrast that we relate to negative charge 
that is localized at the atomic scale is shifted by about 
one lattice constant along the [001] direction. We sug¬ 
gest that this negative charge is the due to the electrons 
that are captured in the broken bonds and this charge 
is shifted away from the Fe atom. Taking into account 
the symmetry of the bright contrast features that are 
related to wavefunction that has been calculated for var¬ 
ious geometries we can conclude that the broken bond is 
aligned along the [001] direction. This corresponds with 
the geometries for the wavefunction calculations shown 
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in Fig gb) or (c). It is quite likely that the presence 
of the s urfac e favors a broken bond along this symmetry 
direct iorP 17 l. 


IV. CONCLUSIONS 

The results of our ab initio calculations show that Fe 
has a second metastable configuration in bulk GaAs with 
a formation energy of 1.2 eV with a configuration similar 7 
to the DX~ center observed in the bistable behavior of 
donors in Al^Gai-^As alloys for x > 0.22. This second 
configuration of Fe differs from the DX~ center by its 
larger displacement from the Ga substitutional position 
to an interstitial position of six-fold coordination sym¬ 
metry. This second metastable configuration is accom¬ 
panied by an associated change of charge state from the 
isoelectronic (Fe 3+ )° state to the (Fe 2+ ) - acceptor state. 
Our results for the case of Si in bulk GaAs also predict 
metastable behavior with an associated change in charge 
state, which is in agreement with the results of previous 
first-principles investigations 1 2 3 4 5 6 7 8 9 ^ of Si in bulk GaAs. The 
tight-binding wavefunction calculations suggest that the 
change in bond configuration required to shift from type I 
to type II Fe does not eliminate the bow tie shape of the 
wavefunction seen experimentally for both types. 

Previous work investigating the bistable behavior of Si 


at the [110] surface of GaAd^has demonstrated that sur¬ 
face effects can be very important in the properties of 
individual dopants, causing very different behavior from 
that seen in bulk. In the case of Si at the surface of 
GaAs, it has been shown in first-principles calculations 6 ^ 
that the formation energy of the DX~ center configu¬ 
ration of Si decreases to become more energetically fa¬ 
vorable than the substitutional configuration if the atom 
resides in the top layer. This work provides evidence that 
a similar behavior occurs for case of Fe when the dopant 
is close to the surface layer and the negative charge state 
is electronically induced by STM. The second metastable 
configuration coordinate of Fe in bulk GaAs along with 
these experimental findings would imply that the bista¬ 
bility behavior of Fe is analogous to the case of Si donors 
in GaAs in which surface effects play a large role in the 
formation of the defect. The bistability of Fe at the sur¬ 
face of GaAs implies the possibility to electronically ma¬ 
nipulate both the charge state and lattice position of a 
TM dopant using STM. 
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